NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 

No. 1859 

A METHOD OF CALCULATING A STABILITY BOUNDARY THAT DEFINES 

A REGION OF SATISFACTORY PERIOD-DAMPING RELATIONSHIP 

OF THE OSCILLATORY MODE OF MOTION 

By Leonard Sternfield and Ordway B. Gates, Jr. 

Langley Aeronautical Laboratory 
Langley Air Force Base, Va. 



Washington 
April. 1949 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE NO. 1859 


A METHOD OF CALCULATING A STABILITY BOUNDARY THAT DEFINES 
A REGION OF SATISFACTORY PERIOD -DAMPING RELATIONSHIP 
OF THE OSCILLATORY MODE OF MOTION 
By Leonard Sternfield and Ordway B. Gates, Jr. 

SUMMARY 


A method has teen derived by which a boundary can be obtained that 
defines a region in which there exists a satisfactory relationship 
between the period and damping of the lateral oscillatory mode of motion, 
according to any given criterion for this relationship. In addition, 
a method is discussed by which curves representing a constant rate of 
spiral divergence may be constructed. 

The methods as presented are applicable to both lateral -stability 
and longitudinal -stability analyses. 


INTRODUCTION 


In lateral-stability analyses it is the usual practice to calculate 
a neutral-oscillatory-stability boundary that is plotted as a function 
of the directional -stability derivative C^ and the effective -dihedral 


derivative C, . This boundary has a special significance, since for a 

l 3 

particular airplane it indicates the combinations of C and C, 

n 3 *3 

necessary for oscillatory stability. Although this boundary is definitely 
an aid in a stability analysis, its value is somewhat restricted in that 
it affords no information about the variation of the period-damping 
relationship throughout the stable region. That is, that region of 
the C n ,Cj plane in which a given criterion for satisfactory relationship 


of the period and damping of the oscillatory mode will be satisfied 
cannot be determined from this boundary. This information can be obtained 
by using the methods described in reference 1 to calculate curves of 
constant period and constant damping, but because of the complexity of 
the calculations involved, a more direct approach to the problem is 
desirable. In the present paper a method is derived by which it is possible 
to obtain directly a boundary, plotted as a function of C n and , " 
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that will define a region in which there exists a satisfactory relation- 
ship between the period and damping of the oscillatory mode of motion, 
according to a prescribed criterion for this relationship.. 

For the combinations of C and C, necessary to satisfy various 

n P 

criterions for the period and damping of th 9 oscillatory mode, undesirable 
spiral instability may be present. In view of this fact a method is also 
presented by which curves representing a constant rate of spiral divergence 
may be constructed. 
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SYMBOIS AND COEFFICIENTS 

I 

angle of roll, radians 
angle of yaw, ■ radians 
angle of sideslip, radians (v/v) 

sideslip velocity along the Y-axis, feet per second 
airspeed, feet per second 
mass density of air, slugs per cubic foot 
dynamic pressure, pounds per square foot 
wing span, feet 
wing area, square feet 
weight of airplane, pounds 
mass of airplane, slugs (W/g) 
acceleration due to gravity, feet per second per second 
relative -density factor (m/pSb) 

inclination of principal longitudinal axis of airplane with respect 
to flight path, positive when principal axis is above flight 
path at the nose, degrees 

angle of flight path to horizontal axis, positive in a climb, degrees 
radius of gyration in roll about principal longitudinal axis, feet 
radius of gyration in yaw about principal vertical axis, feet 
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K- 


-xz 


'n 


^P 


n v 


■'n T . 


nondimens ional radius of gyration in roll about principal 
longitudinal axis ^k x ^/bj 

nondimenaional radius of gyration in yaw about principal 
vertical axis ^k^^/b^ 

nondimenaional radi us' of gyration in roll about longitudinal 
stability axis K x 2 cos 2 ti + 2 sin 2 r) ^ 

nondimenaional radius of gyration in yaw about vertical 
stability axis 2 cos £i T] + K x ^ 2 sin 2 T| ^ 

nondimenaional product-of -inertia parameter 


« 


r 0 o\ 

K z - K x c ) sin Tj cos r\ 


trim lift coefficient ( ff~ coa 7 \ 

\ & ) 

rolling -moment coefficient ^ — 0 11 Ing^moment^ 

yawing -moment coefficient 

lateral -force coefficient ( ^ a , tera, l P° rce ) 

\ <1 S J 

effective -dihedral derivative, rate of change of rolling -moment 
coefficient with angle of sideslip, per radian (dC^dp) 

directional-stability derivative, rate of change of yawing -moment 
coefficient with angle of sideslip, per radian (dC^dp) 

lateral-force derivative, rate of change of lateral-force 
coefficient with angle of sideslip, per radian (dC^/dp) 

damping -in -yaw derivative, rate of change of yawing-moment 
coefficient with yawing -angular -velocity factor, per 


radian 


(*V 


t<. rb y 
2 V 


rate of change of yawing -moment coefficient with rolling -angular - 
velocity factor, per radian 

damping-in-roll derivative, rate of change of rolling -moment 
coefficient with rolling-angular-velocity factor, per 


radian 


(*!/*&) 
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3 b 

Db 

B, 


rate of change of rolling -moment coefficient with yawing -angular 


velocity factor, per radian 


rate of change of lateral-force coefficient with rolling -angular 
velocity factor, per radian 


(**/*&) 


rate of change of lateral-force coefficient with yawing -angular - 


velocity factor, per radian 


rate of change of rolling -moment coefficient with angle of roll, 
per radian ( j<Xj> ) 


time, seconds 


nondimens ional time parameter "based on span (Vt/b) 

0 


differential operator/ — — ) 

\**b) 


Eouth's discriminant 


complex root of stability equation 


AX 


B = \j a 2 + cu 2 

9 = tan" 1 ^ 
a 


+ Bx3 + CX 2 + DX + E = 0 ( 


= a + icjD = Re 


ie ) 


n>r n 


factors used in calculation of boundary of satisfactory period- 
damping. relationship (e n . = - -- , 0 n = — — for 0 < n < M 

y E n E n_1 “ ~ ) 


P 

T l/2 


spiral stability root of stability equation 
AX 1 *' + BX3 + CX 2 + DX + E = 0 
period of oscillation, seconds 

time for amplitude of oscillatory or spiral mode to decrease to 
one -half its original value 

time for amplitude of oscillatory or spiral mode to increase to 
double its original value 


A,B,C,D,E coefficients of lateral -stability equation 
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EQUATIONS OF MOTION 


The nondimensional linearized equations of motion, referred to the 
stability axes, used to calculate the spiral- stability and oscillatory- 
stability boundaries for any flight condition, are: 

Rolling 

+ KxzBb 2 ^)= C 2p p + |cz r Db* + |c Zp Db0 

Yawing 

2m,(K z 2 Db 2 t + KxzIvV) = 0 np S + ^Db* - 
Sideslipping 

SMb^P + V) = + + C I$ + £°Y r + ( C L tan 7^ 

When is substituted for 0, for ^ > and. P 0 e^ 3 ^ for 

in the equations written in determinant form, \ must be a root of the 
stability equation 

AX, 4 + BX 3 + CX 2 + DX + E = 0 (l) 

where 

A = 8u b 3 (K x 2 K z 2 - K XZ 2 ) . 

V 

B = W(s%V% + * *2?% - 

c = l*b( K X 2 c n r c Y|3 + ^b K X 2c np + %, 2 Cl p c Yp + f c n r °Z p * K XZ C l r C Yp 
‘ ^b^IzClp " C n p K XZ C Yp " | C np C 2 r + K XZ C n p C Y p ' K Z 2c Y p C Z p 


- Kx 2 C Yr Cnp + KxzCY r c Zp 
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D = -^ c n r c Z p c Y p ' ^b c 2 p c np + t°np C Z r c Yp + ^b c n p c 2p + 2 ^b c L K XZ c np 
- : 2^t c L K Z 2c Zp - 2u b K x 2c np C L tan 7 + 2 n b K XZ C z p C L tan 7 

+ 1%%°^ - 1%%°^ - 5°I r C n p C Y p + tVV Y P 

E - - %%) + h ^ Y (%% ' %%) 

The damping and. period of the lateral oscillation in seconds are given 
by the equations . 



-0.693 b 
a V 

0-693 b 
a V 

6-28 b 
o> V 


where a and co are the real and imaginary parts of the complex root 
of stability equation (l) . The damping of the spiral mode of motion in 
seconds is given by the equations 


T l /2 


-0.693 b 
X 2 -V 


_ 0 .693 b 

2 " Xg y 


ANALYSIS 


Method, for constructing a boundary that defines a region of satisfactory 
period-damping relationship .- For the derivation of a method by which the 
period- damping relationship of the lateral oscillatory mode may he expressed 
as a function of C n p and C}p, the following assumptions are made: 

1. The stability equation is given in the conventional 

form AX 1 * + BX.3 + CX 2 + DX + E = 0 
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2. The coefficients A, B, C, D, and E are functions of the 
stability derivatives, all of which have prescribed values with the 
exception of C n ^ and C^. 

3- The stability equation has for one of its roots 

X-L = a + ho = Re*® 


where 


0 = tan-1 S 
a 


and 


R = \ia 2 + u£ 


This root Xq represents an oscillatory mode of motion, the period of 
which is given by P = ~ y and the time required for the amplitude of 
the oscillation to d a m p to one-half its original value is given 

^l/2 = y* Substituting Re*® for X in the stability equation, 

and noting that e in ® = cos n0 + i sin n 6 } results in the expression 

AR^(cos 40 + i sin 40) + Br3(cos 30 + i sin 30) 


+ CR2(cos 20 + i sin 20) + DR (cos 0 + i sin 0) + E = 0 (2) 

Equating to zero real and imaginary parts in equation (2) gives the 
equations 

AR^cos 40 + BR 3 cos 3® + Ciscos 20 + DR cos 0 + E = 0 

^ ( 3 ) 

AR^sin 40 + BR 3 sin 3© + Ci^sin 20 + DR sin 0 = 0 

Multiplying the first equation of equations (3) by cos 40, the second 
by sin 40, and adding the two resulting equations give the expression 

Ae 0 + Be^ + C€ 2 + 1*3 + Ee4 = 0 (4) 


where 


6 


n = 


cos n0 

Dp 


Similarly, multiplying the first equation of equations (3) by sin 40, the 
second by -cos 40, and adding the two resulting equations give the 
expression 


A0 O ^ B^-^ 1^3 B^4 = 0 


(5) 
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where 



sin n0 
E n-i 


Since 0 and R are functions of the terms a and cd of the root X-^, 
the values of the factors e n of equation (4) and (f> n of equation (5) 
will he fixed for specified conditions of period and damping. Therefore, 
solving equations (4) and (5) simultaneously results in the 

combination that satisfies any desired period and damping relationship 
of the oscillatory mode. 


If a suitable criterion exists for the damping of the lateral oscillatic 
expressed as a function of the period, for example. 


T l/2 = f ( p ) 


this criterion will be expressible in terms of a and cd by an equation 
of the form 


-0.693 b =f /2n b 

a Y l cd V 


Thus, for any given period, there exists a corresponding value of the 
damping term a necessary to satisfy the criterion. The terms e n 
of equation (4) and terms $n of equation (5) are expr essed in terms 

of R and 0, which are functions of a and cd (r = \ja^~+ ofi 

and 0 = tan'-*- — therefore, curves showing the variation of e n 

and (j> n with the period of the oscillation that will satisfy the given 
criterion can be constructed. For any value of the period P the 
values of ^ and may be taken from these curves and substituted 

in equations (4) and ( 5 ). The corresponding C n ,Cj values may then 

P 3 

be determined from a simultaneous solution of these equations. 


It is important to note that the method described is applicable 
to any criterion for satisfactory damping of the oscillatory mode 
expressed as a function of the period. If the criterion is expressed 
in terms of the number of cycles required to damp to half amplitude, 
that is, = f(P), "the preceding method can still be applied 


since C 


1/2 




General procedure for use of the method . - In the analysis of a 


given airplane the coefficients A, B, C, D, and E of equation (l) are 


evaluated in .terms of the known or prescribed values of all the 
stability derivatives except two /usually C n ^ and which are left in 

symbolic form. A point on a curve y plotted as a function of the two 
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derivatives chosen to he variables, defining the region of satisfactory- 
period- damping relationship is obtained in the following manner: For an 

arbitrary value of the period P, values of e n and $n are taken from 
previously constructed curves that satisfy the desired criterion for 
period- damping relationship of the oscillation. These values are substi- 
tuted in equations (4) and ( 5 ), together with the expressions for the 
coefficients A, B, C, D, and E. The two equations are then solved 
simultaneously for the variable derivatives . By repeating the process 
for successive values of P the desired boundary will be obtained. In 
order that the significance of -the boundary thus obtained be clearly 
understood it Is necessary to know on which side of this boundary there 
exists a satisfactory period- damping relationship according to the pre- 
scribed criterion. A method for obtaining this information is discussed 
in appendix A . 

Application of the method presented to longitudinal- stability analyses .- 
It should be pointed out that the methods presented herein are not limited 
to an analysis of lateral stability but can be readily employed in the 
longitudinal-stability analysis as well. If a suitable criterion exists 
for the period- damping relationship of the longitudinal oscillatory mode 
of motion, curves of 0n and € n plotted against P can be constructed 
and points on a boundary defining the region in which this satisfactory 
relationship exists between the period and damping of the longitudinal 
oscillation can be calculated. The coefficients A, B, C, D. and E to 
be used in the simultaneous solution of equations (4) and ( 5 ) for the 
longitudinal case must, of course, be the coefficients of the longitudinal- 
stability biquadratic. The variables of interest become, instead of C n p 
and Cjp, any two appropriate parameters for longitudinal stability. 

Method for constructing curves of constant rate of spiral divergence .- 
For the combinations of C n p and Cjp necessary to satisfy various cri- 

terions for the period and damping of the oscillatory mode, undesirable 
spiral instability may be present. In view of this fact it is desirable 
to know something about the divergence gradient throughout the C n p,C ^ 

plane . 

The usual method of obtaining this Information involves the solution 
of a series of stability biquadratics, often a rather laborious procedure. 

The following method is presented which gives the desired information with 
a minimum of calculations • 

Assume that the lateral- stability biquadratic (equation (l)) has as 
one of its roots X = X 2 • This substitution for X in equation (l) 
results in the expression 

AX^ + BX 2 3 + CXg 2 + Big + E = 0 (6) 

where the coefficients A, B, C, D, and E are functions of C n and Cjp 
as before. If X 2 = 0, equation (6) Is reduced to the equation P E = 0, 
which defines the boundary of neutral spiral stability. However, if X 2 is 
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assigned a finite value representing a given rate of divergence, equa- 
tion (6) becomes a definite function of and Cjp* The solution of this 


resulting equation for a sequence of values of ^np gives the corresponding 

values of Cj on a curve representing a constant rate of divergence of 

0 

the spiral mode of motion. These curves can be constructed for as many 
values of as may be necessary for an adequate analysis of the 

spiral instability of a particular airplane. 


ILLUSTRATIVE EXAMPLES AND DISCUSSION 


In order to meet the present Navy -Air Force flying-qualities requirement 
(references 2 and 3) the period and damping of the lateral oscillation 
must satisfy the following criterion: 


V = 1 - 5 

(0 g p 1 2) 

T l/2 = 2 ' 5P ' 3 ' 5 

(2 < P < 00 ) 


(7) 


This criterion is illustrated in figure 1. 


As an example of the methods presented, calculations were made for the 
hypothetical high-speed airplane described in table I by using equation (7) 
as the criterion for satisfactory damping of the lateral oscillation. 

e n 

The values of the factors and =- for a sequence of values 

(v/b) n ' (V/b) 0 " 1 

of the period P are presented in table II and the variations of these 
factors with P are shown in figure 2. Note that the ordinates in figure 2 

g „ 0 

are and rather than e and 0 in order that the 

(v/b) n (v/b) n_1 n n 


curves in figure 2 will apply for any airplane at any flight condition if 

the criterion described in equation (7) is used. For a number of values 

of P, the values of e and 0 were determined from these curves and 

n r n 

were then substituted in equations (4) and (5), said the corresponding 
7 values were obtained from a simultaneous solution for these 

•* Q 


V° 


equations. 


The boundary, plotted as a function of C and C, , defining the 

n 0 

region that satisfies the period- damping relationship given by equation (7) 
is presented in figure 3» For purpose of comparison the neutral-oscillatory- 
stability boundary (E^ = 0) for this airplane is also plotted in figure 3. 
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From inspection of the R^ = 0 boundary this airplane is seen to he 

oscillatorily stable for almost any combination of C_ and C 7 shown 

n p *3 

in the figure. This boundary* however* gives no indication of the degree 
of stability present at any point in the quadrant. The curve defining 
the Wavy-Air Force criterion for satisfactory damping of the lateral 
oscillation clearly shows that the criterion is satisfied in only a 
relatively small part of the stable region of the C *C 7 plane. 

n p i p 

Figure therefore* clearly illustrates that the stability of a given 
airplane cannot be accurately evaluated by merely constructing the R^ =. 0 

boundary and noting on which side of this boundary stability exists. 

Additional calculations were made in order to obtain curves of 
constant rate of spiral divergence for this airplane*, and the results are 
presented in figure 4. Should a criterion for the rate of spiral diver- 
gence be established, it will then.be sufficient in a given spiral- 
stability analysis, to plot only the curve representing that rate of 
divergence . 

The roots of the lateral- stability biquadratic usually indicate the 
presence of one periodic and two aperiodic modes of motion. If this be 
true* only one curve defining the region of satisfactory period-damping 
relationship of the periodic mode would exist.. If* however* the roots 
indicate the presence of two periodic modes of motion* two or more branches 
of this curve may exist. 


Thus* reference 4 shows that for an airplane equipped with an 
automatic pilot that applies aileron control proportional to the dis- 
placement in roll, the roots of the stability equation will represent two 
oscillatory modes of motion. Calculations were made for this configuration 
by using the stability derivatives and mass characteristics presented in 
table I of reference 4. Points on the boundary of satisfactory period-- 
damping relationship were calculated for C 7 = -0.10 by using as 

ip 

variables the directional- stability derivative C and the deriva- 

n P 

tive C* due to the automatic pilot. The equations and coefficients 

0 

presented in the section entitled "Equations of Motion" were modified to 


take into account the derivative C- 


^0 


in the same manner as was done in 


reference 4. 


The results of these calculations are presented in figure 5 ° For a 
number of points in the C-q^* plane* roots of the stability equation 


lp 
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were calculated in an effort to understand more clearly the significance 
of each of the tranches of the boundary presented in this figure. The 
results of these calculations are given in table III. At point A the 
roots represent two periodic modes, . neither satisfying the Navy-Air Force 
criterion. Upon passing through the branch of the boundary in the upper 
quadrant to point B, the period- damping relationship of one of the 
periodic modes becomes satisfactory while the other remains unsatisfactory. 
At point C the relationship of the unsatisfactory modes has improved some- 
what and the other mode remains satisfactory. At point D the previously 
satisfactory periodic mode has become two aperiodic modes, one of which 
is unstable. The period- damping relationship of the, periodic mode is now 
satisfactory. At point E the periodic mode has again become unsatisfactory 
and the two aperiodic modes are stable. Thus, the only part of 
the C n ,C^ plane that completely satisfies the Navy- Air Force criterion 


is that area on the unhatched side of the branch of the boundary appearing 
in the lower quadrant. It should be pointed out that although the 
oscillatory mode satisfies the Navy-Air Force criterion within this area, 
the instability of one of the aperiodic modes is such that the region 
may be of little practical value. 


CONCLUSIONS 


The following conclusions were made from a theoretical investigation 
to develop a method of calculating a stability boundary that defines a 
region of satisfactory period-damping relationship of the oscillatory 
mode of motion: 

1. Through use of the methods presented a boundary can be obtained 
that defines a region in which there exists a satisfactory relationship 
between the period and damping of the lateral oscillatory mode of motion, 
according to any given criterion for this relationship. 

2. A method is also presented by which curves representing a constant 
rate of spiral divergence may be constructed. 

3. These methods, which are developed in detail herein for the 
analysis of lateral stability, are adaptable as well to the analysis of 
longitudinal stability. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., February 23, 1$&9 
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APPENDIX A 

A METHOD FOR DETERMINING ON WHICH SIDE OF THE BOUNDARY 
OF SATISFACTORY PERIOD- DAMPING RELATIONSHIP 
A GIVEN CRITERION IS SATISFIED 


If the stability biquadratic has for one of its roots X^ = a + Ixu, 
it can be shown from reference 1 that the following parametric equations 
are satisfied 


f-L - nf 3 = 0 (Al) 

f - [if 2 + P 2f '4 =0 (A2 ) 

where ' " ■ • • 

[ 1 = 0 ? 

f --= Aa^ + Ba 3 + Ca 2 + Da + E 
f^ = 4Aa 3 + 3 Ba 2 + 2Ca + D 


f 2 = 6 Aa 2 + 3Ba + C 
f 3 = 4Aa + B 
f 4 = A 

The coefficients A, B, C, D, and E appearing in the expressions for f, 
f]_, f 2 , f 3 , and f^ are those of the characteristic stability equation. 
These coefficients are functions of the stability derivatives, two of 
which are assumed to be variable . 

If the parameter n Is eliminated between equations (Al) and (A2), 
the expression obtained is 


(f 2 f 3 - fif4)fi - f 3 2 f - 0 


(A3) 


Thi 3 equation represents a curve of constant damping a, and the frequency 
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of the oscillation at any point on the curve can "be calculated, by using 
equation (Al) . If a = 0, equation (A3) reduces to the familiar expres- 
sion for the neutral-oscillatory- stability boundary (BC - AD)D - B 2 E = 0. 

Now, for every point on the boundary of satisfactory period and 
damping, plotted as a function of any two arbitrary derivatives x and y, 
a root = a + will exist that exactly satisfies the prescribed 

criterion. Assume that at a particular point on the boundary (xi,y^) the 
root is Xi = ai + hoi- If a = ai and x = xg = x^ + Ax be substi- 
tuted in equation (A3), the value 7 - 72 ) which is located- on a curve 
of constant damping, can be calculated. By substituting aj_, xg, and 72 
for a, x, and y, respectively, in equation (Al) the value of n, and 
hence 0)9, is determined at the point (xg, 72 )' If the root a^ +'im 2 at 
point (xg,y 2 ) satisfies the prescribed criterion, the region in which 
this point is located is the satisfactory region with respect to the 
boundary defining the criterion. 

An imaginary value of aug at point (xg,y2) indicates that no com- 
plex root with the real part equal to a^ exists at this point, since cog 
must be real if the root ai + icu2 la to represent an oscillation. 

That is,- a^_ + lug will represent two real roots and before the 
point (xg, y 2 ) can be established as satisfactory or unsatisfactory, the 
other two roots of the stability equation must be determined. However, 
if Ax ia chosen small enough, ag will be real. 

In general, the satisfactoiy region may also be readily identified 
if roots of the stability biquadratic are calculated at several points 
of interest. A method of evaluating the roots of a quartic equation is 
presented in appendix B . 
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APPENDIX B 


METHODS FOE EVALUATING THE BOOTS OF A QCJAETIC EQUATION 


Various methods exist for evaluating the roots of high-order poly- 
nomials (references 5 and 6) . Many of these methods, although highly 
accurate, become rather laborious in actual application. 

A method developed by Lin (reference 5) and independently by Doris 
Cohen of the NACA (unpublished) affords a means of evaluating roots of 
high- order equations with a minimum of computation. This method is 
basically one of synthetic division and under certain conditions, which 
are discussed subsequently, the desired roots can be obtained very 
rapidly • 

In order to illustrate the method consider the quartic equation 


\ k + BA 3 + CX 2 + EX + E = 0 

As a first approximation, assume that the equation has the quadratic 

factor A 2 + JX + k where J = — and k = — . Division may then be 

C C 

perfoimed as follows: 


A 2 


+ JX + k 


X 2 + (B - jjX + C - k - j(B - j) 
+ BA 3 -* CX 2 + DX + E 


(C - k - J(B - JI|X 2 + 
[c - k - J(B - J)]a 2 + 


[5 - k(B - jj| X + E 

j jc - k - J(B - jj X + kg - k - J(B - jj ] 
Bemainder 


If the remainder is equal to zero or negligible, the two quadratics 
obtained are good approximations to the factors of the quartic equation. 
If the remainder is not, negligible, the procedure is repeated by using 
as the second approximation the factor X 2 + j'X + k' where 
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and 


k ' C - k - J (b - j) 

This operation is continued until the remainder is negligible . The results 
usually converge rapidly if the frequencies of the two modes are of differ- 
ent orders of. magnitude or if the- quartic equation has small real roots . 

The usefulness of the method is greatly reduced in cases where the results 
converge slowly. 


A substitute method has been derived by the authors to be used as a 
means of obtaining the roots of the quartic equation if the. method of 
reference 5 converges slowly . 


It can be seen from the division performed previously that 
if X 2 + JX, + k is a factor of the quartic equation, two parametric equa- 
tions must be satisfied, namely 


— \ 

D - k(B - J) = j[c - k - J (B - j]] 


(Bl) 


E = k[c - k - j (B - jJ 

The method consists in solving these two equations simultaneously for J 
and k. A convenient procedure is as follows: 


Equations (Bl) are rearranged in the form 

k _ J3 - jgBt JC - D = Q 1 
2J-B 

( 

k 2 - (j 2 - JB + C)k + E = 0 


(B2) 


The first of equations (B2) is readily solved for k for any given 
values of J. A series of corresponding values of j and k are computed 
in this way. For each corresponding j and k the value of the left- 
hand side of the second equation of equations (B2) is evaluated. This 
value is plotted as ordinate against j as abscissa. The intersection 
of the curve with the j-axis evidently determines the solution for J . 

The solution for k is obtained upon substitution of this value of j 
back into the upper equation. The values of j and k thus obtained 
may be substituted back into the. lower equation to check the solution. 


The methods presented are readily adaptable to equations of any even 
order. In the case where an uneven- order equation should appear, the 
logical procedure is to obtain a real root by synthetic division and there 
by reduce the equation to an even order. 
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TABLE I 

STABILITY DERIVATIVES AND MASS CHARACTERISTICS 
OF HYPOTHETICAL AIRPLANE 


W/S, lb /ft 2 

b,.ft . . . . 

p, slugs/ft^ . • . '. . . 
V, ft/sec 

7, deg . . . . . . • ... . . . 

C L * * • ’ * * % . * *• • • 

^b 

k Xo,ft 

k z 0 / ft; ... . 
T l, deg. 

C Zp , per radian 

C^’/ pfer. radian . . . ' . . 
C np , per radian . . . . . 

c n r > P er radian 

Cy p , per radian 

c Y r > P er rad -ian 

Cy^, Per radian 

%(fuselage)’ radi “ 
C n p, P er radian 

Cn p(tail ) ’ per radian • * 
C^p, per radian 



. . . 40 
. .. . 20 

; 0.0004 

... 733 
... 0 

. . 0.372 

, . . 155 
* • 2.02 
; 9.64 
... 2 

. . - 0.197 

0.0929 
-0 .00732 


-1 


-1 


47 


33 


'Hp (tail) 

. . . 0 

. . . 0 

'' n P(tail) 

. . -0.25 


Variable 

° n p ‘ ° n p( fuselage) 
Variable ‘ 
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TABLE III 

NATURE OF ROOTS OF STABILITY EQUATION 
™ C np ,C PLANE 


Point 

(see fig. 5) 

T l/2 

T 

2 

p 

A 

IO.65 


4.95 


3*08 

* “ — “ 

1.19 

B 

9-95 

— 

7.05 


3*16 

. 

1.53 

C 

7-72 

— 

30.2 


3-50 

- - “ - 

2.16 

D 

1.02 

1.65 

Aperiodic 


4.20 



3-15 

E 

.793 

2 *5P_ 


— 

Aperiodic 



2.39 

7-93 
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Figure !•— Navy-Air Force criterion for period— damping relationship of the oscillatory mode of motion. 


Figure 2 Variation of the computation factors 


with 
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S — and. 

(V/b) n 

for the Navy-Air Force damping-period, criterion. 


and — — - are presented in table II. 


. Fn 

(v/b ) n— 1 

The values 







Figure 3 .— Boundary for satisfactory relation "between damping and period according to Navy-Air Force 
criterion. R-, = 0 boundary shown for comparison. Hypothetical airplane. 






Figure 5.— Boundary defining region of satisfactory period— damping 
relationship according to Navy-Air Force criterion. Airplane with 
automatic— pilot coupling of aileron to angle of hank. 





